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Using a previous extraction of the quark Boer-Mulders distributions from semiinclusive deep 
inelastic scattering data, we fit the unpolarized Drell-Yan data on the cos 2<j> asymmetry, determining 
the antiquark Boer-Mulders distributions. A good agreement with the data is found in the region 
of low qr, where the transverse- momentum factorization approach applies. 
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Ifl \ I. INTRODUCTION 

r-> : 

One of the most relevant results of high-energy spin physics in the last decade has been the discovery of many 
interesting correlations between the transverse momentum and the transverse spin of quarks (for reviews, see Ref. [l|, 
. 0]). A surprising consequence of these correlations is that there may exist non trivial spin effects in unpolarized hard 
O i 1 processes, generated by a leading- twist, chiral-odd, transverse-momentum dependent distribution function, the so 
called Boer-Mulders function h^(x,\t?p) [|[, which represents a transverse-polarization asymmetry of quarks inside 
qj- an unpolarized hadron. The origin of was clarified in Refs. [3-0] and the first calculation in a realistic quark- 
r~| \ diquark model was reported by Goldstein and Gamberg ||. In 1999 Boer H s ug g ested that could explain the 
1-1 . large cos 2(f) asymmetries observed in unpolarized irN Drell-Yan production [1CH12| , which were not understood in 
terms of purely perturbative QCD effects [HI EH H] ■ This finding was confirmed by more refined model calculations 
[Ti| . A comparable or even larger asymmetry is predicted for pp D rell-Yan production [l6l - [2lj , a process to be studied 
£f~) • in the next years at the GSI High-Energy Storage Ring [22l. l23l| . While nN and pp probe valence distributions, pp 
£SJ , and pD Drell-Yan reactions are sensitive to the sea distributions and therefore the corresponding cos 2(f> asymmetries 
■"sj" ■ are expected to be smaller. This is indeed what the E866/NuSea experiment found: the cos 2<f> dependence observed 
C^) | in pD dimuon production is of the order of few percent [24 . [25| . 

. A cos2(/) asymmetry also occurs in unpolarized semiinclusive deep inelastic scattering (SIDIS), where it has been 
measured in the low transverse- momentum region by HERMES [H|, COMPASS [27| and CLAS [H|. In SIDIS the 
Boer-Mulders distribution couples to a chiral-odd fragmentation function, the Collins function H ^ [29|, which de- 
scribes the fragmentation of transversely polarized quarks into polarized hadrons. Recently, we presented a systematic 
phenomenological analysis of the various contributions to the cos 2(f) asymmetries in unpolarized SIDIS [30j , and of the 
preliminary HERMES and COMPASS results [3l|. In the kinematics of these experiments the perturbative term is 
found to be negligible, whereas the order-fc^/Q 2 contribution from non-collinear kinematics (the so-called Cahn effect 
[32|,[33|) is quite large. The Boer-Mulders effect is also sizable and generates a negative (positive) asymmetry for n + 
(tt~), due to the expected negative sign of the u distribution. Combining the Cahn contribution (which is positive 
and roughly the same for 7r + and tt~ ) with the Boer-Mulders contribution, we predicted a n~ asymmetry larger than 
the 7r + asymmetry [30| . Our results turned out to be in fair agreement with the first SIDIS measurements of the 
cos 2(f> asymmetry, as shown in our most recent study [3l| , where we attempted an extraction of the Boer-Mulders 
distributions from the HERMES and COMPASS preliminary data. Since the present statistics is not sufficient to 
allow a complete determination of , the strategy of Ref. [3l| was to use the same functional form as the Sivers 



function obtained from SIDIS data in Ref. [34J and parametrize the normalization coefficient. Results compatible 



with impact-parameter expectations [351 ] combined with lattice findings |36j . and with model calculations 37[ were 
found. However, the SIDIS data leave the Boer-Mulders sea totally unconstrained. Thus, in order to get the antiquark 
distributions hj~ one needs an extra source of information. 

The purpose of the present paper is indeed to extract some information on the antiquark Boer-Mulders distributions 
from the Drell-Yan pp and pD data on the cos2</> asymmetry [U, [25j]. We perform a fit to these data in the region 
of low qt, where the transverse- momentum factorization is expected to be valid [38[ and the perturbative effects are 
small. Since in Drell-Yan processes the Cahn contribution to the cos 2<fi asymmetry is known to be negligible [39[ , an 
analysis in terms of the Boer-Mulders effect alone is possible. 
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II. THE cos 2cf> ASYMMETRY IN DRELL-YAN PROCESSES 

The angular differential cross section for the unpolarized Drell-Yan process is usually parametrized as 
1 da DY 3 1 



, 1 + A cos 2 9 + u sin 29 cos 6 + - sin 2 9 cos 26) . (1) 
4ttA + 3V r ^ 2 / w 

where 6 and 6 are > respectively, the polar angle and the azimuthal angle of dilcptons in a dilcpton center of mass 
frame. In particular, we adopt the Collins-Soper frame (40| . where 9 is the angle between the dilepton axis and the 
bisector of Pi and — P2 (the momenta of the colliding hadrons), and 6 is the angle between the lepton and hadron 
planes. We denote by qx = \qt\ the transverse momentum of the lepton pair (or, equivalently, of the virtual photon). 
The v parameter in JT]) is the cos 26 asymmetry we are interested in. 

A non-collinear factorization theorem for Drell-Yan process has been proven by Ji, Ma and Yuan [38| for qr -C Q 
(where Q is the invariant mass of the lepton pair). At order o^, the 0- independent term of the unpolarized Drell-Yan 
cross section is 



dfl dx\ dx2 d 2 qx 12 Q 2 



(l+cos 2 0) J d 2 k 1T d 2 k 2T 6 2 (k 1T +k 2T -q T )[f^x 1 ,k 2 1T )^(x 2 ,k 2 2T ) + (l ^ 2)}, . (2) 



Here /i(cc,k 2 T ) is the unintcgrated quark number density. 

The Boer-Mulders contribution to the unpolarized cross-section reads @ 



da 



DY 



dft dx\ dx 2 d 2 qx 



a 



2 

cm _:„2 



cos 2<p 
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1" 9 ^2 e l J d 2 ki T d 2 k 2 T 5 2 {ki T + k 2T - q T ) 



(2 h • k iT h • k 2 r - ki T • k 2T) 



m 



[hf a {x 1 ,ki T )hf a (x 2 ,k 2 T )cos26 + (1«2)], (3) 



(4) 



with h = qx/qr- 

From Eqs. © and © we get the following expression for the coefficient v (setting A = 1, fj, = 0) : 

with the following notations: 

/«] - J d 2 k 1T d 2 k 2T S 2 (k 1T + k 2T - q T ) x f?(x u fc? T )/f (x 2 , k 2 T ) 

= Jdk 1T k 1T ^dx/ 1 a (^i,^T)/f(^2,|qT-k 1T | 2 ), (5) 

n[hi a ,hi a ] = [ d 2 k 1T d 2 k 2T 5 2 (k 1T + k 2T - q T ) (2h ' klTh ' k2T 2 ~ klT ' k2T) h^( Xl , k 2 1T )h^(x 2 , k\ T ) 

f ai 1 f 27T a k 2 T + q T k 1T cos X ~ 2 k 2 T cos 2 x ,j_ af , 2 uia, 1 , \2\ fR \ 
= dk 1T k 1T / d X — - 2 - hf a (x 1 ,kf T )hf a {x 2 ,\q T - k 1T | 2 ) , (6) 

J JO m N 

where \ is the angle between qr and ki^ and we omitted for simplicity the (1 O 2) terms. The asymmetry v in 
Eq. ([4J depends on the kinematic variables X\, x 2 , Q and qx- 

In general, there is another contribution to v arising from the Cahn effect, that is from purely kinematic transverse- 
momentum corrections to the ordinary parton model formulas. However, as shown in Ref. [39j . the Cahn contribution 
to v is proportional to {{k 2 T ) — (k 2 T )) 2 and hence negligible, or even strictly vanishing when the average transverse 
momenta of quarks (or antiquarks) in the two colliding hadrons are equal (which is what we assume here). Thus, the 
Boer-Mulders effect is the only non-perturbative source of a cos 20 asymmetry up to order q\IQ 2 ■ 
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N u = -18 


N d 


= -45 


a u = 0.73 




= 1.08 


Pu=P d = 3.46 


M 2 


= 0.34 GeV 2 



TABLE I: Parameters of the Boer-Mulders quark distributions [3"t| . 



III. PARAMETRIZATIONS OF DISTRIBUTION AND FRAGMENTATION FUNCTIONS 

Let us consider first of all the /c^-dependent unpolarizcd distribution functions (where kq< stands either for k\T or 
for /c2t)- We assume that these functions have a Gaussian behavior in fcr, 

p—k T / {k T ) 

tf{x,h 2 T ) = tf{x) (7) 



which is supported by lattice studies [41( and by a recent phcnomenological study of SIDIS and DY cross sections 
[39| . The integrated unpolarizcd distribution functions /' are taken from the GRV98 fit 42 1. 



Since the available SIDIS data do not allow a full extraction of the Boer-Mulders function, in Ref. [3l| we assumed 
to be simply proportional to the Sivers function f^,, 

hi a (x,k 2 T ) = X a f 1 A T a (x,k 2 T ), (8) 



with f^f taken from a phenomenological analysis of the Sivers asymmetry [34| and the coefficient A a fitted to the 



SIDIS cos 20 data. Various theoretical arguments (based on the impact-parameter picture |35j , on large- N c arguments 
[43I] . and on model calculations [HI, |45[) suggest that the u and d components of hj~, at variance with should 
have the same sign and in particular be both negative (which means that Xd should be negative). This is indeed what 
we found in Ref. [3l[. Moreover, the impact-parameter approach [35[ combined with lattice results [36| predicts a u 
component of hi larger in magnitude than the corresponding component of f^, and the d components of hi and 
/j^ with approximately the same magnitude (and opposite sign). 

In our SIDIS analysis the quark Boer-Mulders distributions were parametrized as 

hi a {x,k 2 T ) = N a x^(l~x)^ C - k ^^ 2 f1{x,k 2 T ) (9) 

with the parameters a a ,/3 a ,/i borrowed from Ref. [H| and N a fitted to the data. The values of these parameters are 
collected in Table HI Notice that^when used to calculate DY asymmetries, the Boer- Mulders distributions determined 
in SIDIS must be sign-reversed [J, [f| . 

Concerning the Boer-Mulders antiquark distributions, the SIDIS data are quite insensitive to them. Therefore, in 
Ref. [3l| these distributions were not fitted to the data, but just taken to be equal in magnitude to the corresponding 
Sivers distributions. 

The Drell-Yan process, on the contrary, probes the products of quark and antiquark distributions. Thus the Drell- 
Yan measurements of the cos 2<j> asymmetry [24], [25| can at least in principle give information about the antiquark 
sector of the Boer-Mulders function. An analysis of the DY cos 2<f> asymmetry has been already performed by other 
authors 0, [47|. However, when extracting the Boer-Mulders distributions from the present DY data, one should 
keep in mind that the fc^-factorization approach applies to the low- qj r egion only, whereas at large qx the observed 
v values are likely to be explainable in terms of perturbative QCD [48l |49(. Therefore, in our analysis we will only 
consider the low-gr DY data. 

The antiquark Boer-Mulders distributions are parametrized as 

hi a (x,k*)=N a x a *(l-x)e*e- k T^ 2 /?0r,fc|), (10) 

with the same a, j3 and (i parameters as for the Sivers antiquark distributions, and the normalization coefficients N a 
and iVj fitted to the data. 

The final parameters to be considered are the average values of the quark and antiquark transverse momenta in 
the two hadrons, (k 2 T ) and {k\ T ). We take these parameters to be equal to each other, {k 2 -) = (k 2 T ) = (fc| T ), and 
we choose for them two different values: 



Fit 1 : = 0.25 GeV 2 ; Fit 2 : (k 2 ,) = 0.64 GeV 2 . 



(11) 
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FIG. 1: The curves represent the results of our fit I to the v asymmetry of pD DY production. Data are from Ref. [2J 



The smaller value, 0.25 GeV 2 , is the one we used in our analysis of the Boer-Mulders effect in SIDIS [31(, and is taken 
from the study of the Cahn effect of Ref. [50[ • 

The larger value, 0.64 GeV 2 , is the one obtained by D'Alesio and Murgia [5l| in their analysis of pp scattering data 
at y/s ~ 20 GeV. A recent phenomenological study of the transverse-momentum dependence of DY cross sections has 
obtained a very similar value [39| . 



IV. RESULTS 

The E866/NuSea Collaboration presented data on the angular distributions of DY dimuons in pD [24| and pp 
interactions [25[ over the kinematic range: 

4.5<Q<15GcV, 0<q T <4GeV, 0<x F <0.8. 

Concerning q?, as already mentioned, we will only consider the data in the q?> < 1.5 GeV region, where the kx 
factorization applies and the perturbative contribution, roughly proportional to q T /Q 2 , is small. 

Our fit 1 to the E866/NuSea data on v is shown in Fig. [T] for pD and in Fig. [2] for pp. The x 2 P er degree of 
freedom is 1.24. The vertical error bars represent the statistical uncertainties. Notice that the experimental bins in all 
variables (xi, X2 and qr) are quite large: for instance, the three points in the qr plot correspond (in GeV) to (0,0.5), 
(0.5,1.0), (1.0,1.5). In all figures we only showed the central values of the bins. 

Concerning fit 2, the corresponding curves overlap those of fit 1, with the same value of x 2 /d.o.f. The difference 
between the Gaussian widths of the distribution is in fact compensated by the different normalizations of the antiquark 
distributions that we obtain from the two fits. 

The values of the parameters of the antiquark distributions are listed in Table [Til The first k\ moments of the 
antiquark distributions, i.e. 



N 

are plotted in Fig. [3] One may notice the sensible difference between the distributions extracted from the two fits. 
The width of the Gaussian distribution is clearly a crucial parameter, but the scarcity of present data docs not allow 
extracting it from the fit. A combined analysis of SIDIS and DY data, and possibly of various azimuthal asymmetries, 
might help to improve the situation. 

Finally, let us compare the present work with with the results of a recent analysis [25| of the E866/NuSea measure- 



ments. In Ref. 25| both the quark and the antiquark distributions have been extracted from pp and pD DY data. 
Since in the DY cross section hj~ a always couples to h^ a , only the magnitude of the products h^ a hj- a is actually de- 
termined, and the signs of the distributions are undefined. In our analysis the use of both the DY and the SIDIS data 
allows constraining separately the magnitudes and the signs of the Boer-Mulders quark and antiquark distributions. 
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FIG. 2: Our fits to the v asymmetry of pp DY production. Data are from Ref. [25 



a* = a d = 0.79, /So = fa = 3.46, fi 2 = 0.34 GeV 2 


Fit 1 


JVs = 3.6 ± 1.0 


Nj = 1.7 ± 1.4 


(k 2 T ) = 0.25 GeV 2 


Fit 2 


N-a = 6.4 ± 1.7 


N s = 3.0 ± 2.4 


(fer) = 0.64 GeV 2 



TABLE II: Parameters of the Boer-Mulders antiquark distributions. 



V. CONCLUSIONS 

We performed a fit to the DY data on the cos 20 asymmetry, showing that using the Boer- Mulders quark distri- 
butions previously extracted from SIDIS and a new set of antiquark distributions one can obtain a reasonably good 
description of v. However, the resulting Boer-Mulders functions depend rather strongly on the width of the Gaussian 
distributions, that is on the average k\. In order to achieve a more precise determination of h\ for quarks and 
antiquarks a combined fit of SIDIS and DY data, and of cos </> and cos 2<f> asymmetries, must be performed. This work 
is now in progress. 
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FIG. 3: The first k% moments of the antiquark distributions from Fit 1 (solid curves) and Fit 2 (dashed curves). 
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